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Microwave  near-field  quantum  control  of  spin  and  motional  degrees  of  freedom  of  ^^Mg+  ions  can  be  used 
to  generate  two-ion  entanglement,  as  recently  demonstrated  in  Ospelkaus  et  al  [Nature  476,  181  (2011)].  Here, 
we  describe  additional  details  of  the  setup  and  calibration  procedures  for  these  experiments.  We  discuss  the 
design  and  characteristics  of  the  surface-electrode  trap  and  the  microwave  system  and  compare  experimental 
measurements  of  the  microwave  near  fields  with  numerical  simulations.  Additionally,  we  present  a  method  that 
utilizes  oscillating  magnetic-field  gradients  to  detect  micromotion  induced  by  the  ponderomotive  radio-frequency 
potential  in  linear  traps.  Einally,  we  discuss  the  present  limitations  of  microwave-driven  two-ion  entangling  gates 
in  our  system. 
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I.  INTRODUCTION 

Scaling  to  large  numbers  of  quantum  bits  (qubits)  is 
important  for  the  development  of  a  future  quantum  information 
processor.  Several  pathways  toward  this  goal  have  been 
proposed  based  on  using  trapped  atomic  ions  as  qubits  [1-6]. 
One  particular  approach  envisions  processing,  transporting, 
and  sympathetic  cooling  of  ions  inside  a  large  array  of  traps 
[2-4].  Recent  experiments  with  various  trap  architectures  have 
demonstrated  the  key  ingredients  for  scalable  ion  loading, 
trapping,  processing,  and  transport  [7-17]. 

Current  implementations  of  quantum  information  proces¬ 
sors  require  a  reduction  of  gate  errors  to  be  efficiently  scalable 
[18],  necessitating  stringent  qubit  control.  In  most  trapped-ion 
quantum  information  experiments,  control  is  accomplished  via 
laser-based  techniques  [4],  but  other  approaches  are  being 
investigated,  e.g.,  those  based  on  magnetic  fields  [19-21]. 
These  methods  reduce  the  required  laser  overhead  and  may 
enable  a  higher  level  of  control  and  integration  than  laser- 
based  control.  While  the  highest  entangled- state  fidelity  is 
currently  obtained  with  laser-based  operations  [22],  magnetic- 
field-based  control  may  eventually  provide  sufficiently  low 
errors  for  fault-tolerant  operations  [23-27].  For  single-qubit 
operations  performed  with  microwave  magnetic  near  fields, 
errors  per  computational  gate  of  2  x  10“^  and  better  have 
been  reported  [28,29].  Although  the  theoretical  threshold  for 
fault  tolerance  depends  strongly  on  the  assumptions  of  the 
model  under  study,  an  operational  error  of  10“"^  or  below 
may  be  a  reasonable  starting  point  in  a  future  fault-tolerant 
architecture.  Reaching  this  level  of  precision  for  two-qubit 
operations  requires  improvements  on  the  state-of-the-art  errors 
by  approximately  two  orders  of  magnitude. 
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As  with  conventional  computers,  microfabrication  may  be 
advantageous  in  realizing  a  future  quantum  processor,  such 
as  those  based  on  surface-electrode  ion  traps  [30]  (for  a 
review  see  Ref.  [31]).  In  addition  to  reducing  the  overall 
volume  of  a  processor  for  a  given  number  of  qubits,  a 
reduction  in  scale  could  allow  for  tighter  confinement  with 
trap  frequencies  scaling  approximately  as  1/J^  if  all  other 
parameters  are  held  constant.  Here  d  denotes  a  characteristic 
length  scale  of  the  trap,  often  characterized  by  the  smallest 
ion-electrode  distance.  Higher  trap  frequencies  can  enable 
shorter  durations  for  transport,  separation  and  recombination, 
and  two-qubit  operations  mediated  by  ion  motion.  In  the  past 
two  decades,  typical  trap  length  scales  have  been  reduced 
from  centimeters  to  tens  of  micrometers.  This  increases  the 
oscillating  magnetic  near  fields  (scaling  as  1/J),  originating 
from  oscillating  currents  running  in  one  or  multiple  trap 
electrodes.  Fast  and  precise  single-qubit  rotations  (carrier 
transitions)  have  been  demonstrated  in  this  way  [28,29,32]. 
In  addition,  the  corresponding  field  gradients  (proportional  to 
l/d^)  can  become  large  enough  to  couple  the  ions’  spin  with 
their  motional  state  (sideband  transitions)  on  experimentally 
accessible  time  scales  [21,32].  For  trapped  ions  close  to  their 
motional  ground  states,  the  ratio  of  sideband- transition  rates 
to  carrier-transition  rates  driven  by  oscillating  near  fields  and 
their  corresponding  field  gradients  is  on  the  order  of  a^/d, 
where  is  the  harmonic-oscillator  ground- state  extent  [21]. 
For  typical  experimental  parameters,  uq  is  on  the  order  of 
10  nm,  and,  therefore,  a^/d  ranges  between  10“^  and  10“"^  for 
the  smallest  traps  today.  Further  reduction  in  trap  size  would 
increase  this  ratio,  enabling  faster  spin-motional  entanglement. 
For  the  same  reason,  near-field  gradients  are  typically  not 
useful  for  larger  (J  >  100  /xm)  trap  structures. 

The  main  practical  limitation  to  size  reduction  has  been 
“anomalous”  motional  heating,  caused  by  electric-field  noise 
at  the  positions  of  the  ions  in  excess  of  Johnson  noise  from 
the  resistances  of  the  electrodes  and  of  their  filters  [33,34]. 
Anomalous  noise  fields  have  been  phenomenologically  mod¬ 
eled  as  originating  from  potential  fluctuations  of  surface 
patches  or  dipoles  with  typical  dimensions  much  smaller  than 
J,  which  yields  a  noise  spectral  density  that  scales  as  1/ [33] . 
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As  d  is  decreased,  the  resulting  heating  grows  faster  than  trap 
frequencies,  oscillating  near  fields,  and  their  corresponding 
gradients.  Consequently,  this  form  of  heating  appears  to 
prevent  scaling  to  even  smaller  trap  dimensions.  However, 
recent  studies  have  demonstrated  a  substantial  reduction  of 
anomalous  heating  by  cooling  the  electrodes  [35-37]  and  by 
treatment  of  the  electrode  surfaces  [38,39].  Such  methods  may 
enable  the  further  reduction  in  trap  size. 

In  the  experiments  described  here,  we  investigate  building 
blocks  for  oscillating  magnetic  near-field  quantum  control  of 
trapped  ions  [2 1  ] ;  a  similar  approach  is  pursued  in  Ref.  [40] .  To 
test  the  basic  scheme,  we  integrated  three  microwave  current- 
carrying  electrodes  into  a  room-temperature  surface-electrode 
trap  with  J  30  /xm.  Frequencies  of  the  microwave  currents 
are  tuned  to  near  resonance  with  hyperfine  transitions  in  ^^Mg+ 
(nuclear  spin  5/2).  The  resulting  magnetic  field  at  the  position 
of  the  ions  determine  the  carrier  coupling,  while  the  magnetic- 
field  gradient  couples  the  hyperfine  levels  (spin  states)  to  the 
motional  degrees  of  freedom  [21].  We  previously  reported  on 
fast-carrier  transitions  and  generation  of  two-ion  entanglement 
[32];  the  entanglement  was  produced  with  a  two-qubit  gate 
operation  [41-43]  suitable  for  scalable  quantum  information 
processing.  In  this  report,  we  give  a  more  detailed  account  of 
our  experimental  setup  and  the  procedures  for  calibrating  and 
optimizing  the  operation  of  the  system. 

In  Sec.  II,  we  discuss  design,  fabrication,  and  performance 
of  the  surface-electrode  trap  and  describe  the  laser  beams  for 
loading,  state  preparation,  and  detection.  Section  III  describes 
the  microwave  system  that  generates  and  shapes  the  signals 
used  for  coherent  single-qubit  and  two-qubit  operations.  In 
Sec.  IV,  we  give  an  experimental  sequence  for  calibrating 
and  studying  entangling-gate  operations.  The  experimental 
procedure  to  adjust  the  microwave  field  at  the  ion’s  position 
for  driving  sideband  transitions,  a  prerequisite  for  entangling 
gates  [41-43],  is  described  in  Sec.  IV  A.  Section  IV B  details 
a  technique  for  rf  micromotion  detection  by  observing  tran¬ 
sitions  in  microwave  magnetic-field  gradients.  In  Sec.  IV  C, 
we  describe  a  procedure  to  optimize  the  orientation  of  the 
radial-motional  modes  relative  to  the  microwave  field  gradient 
components  and  a  static  homogeneous  quantization  field. 
Finally,  in  Sec.  V  we  discuss  the  factors  currently  limiting 
two-qubit  gate  fidelity  [32]. 


II.  SURFACE-ELECTRODE  TRAP  AND 
EXPERIMENTAL  SETUP 

The  surface-electrode  trap  is  designed  as  an  asymmetric 
“five-wire”  trap  [44,45].  It  consists  of  two  rf,  six  control,  and 
three  microwave  electrodes;  see  Fig.  1  for  a  view  of  the  trap 
assembly  [Figs.  1(a)  and  1(b)]  and  the  central  region  of  the 
trap  electrodes  [Fig.  1(c)]. 

The  trap  chip  is  fabricated  in  the  NIST  Boulder  cleanroom. 
We  choose  aluminum  nitride  (AIN)  as  a  substrate  for  its  high 
thermal  conductivity  of  approximately  170  W/(m  K),  which 
aids  dissipation  of  heat  generated  by  the  microwave  currents  in 
the  trap  electrodes.  Adhesion  (Ti,  10  nm)  and  seed  (Cu,  60  nm) 
layers  are  deposited  on  a  75-mm-diameter  AIN  wafer  by 
electron  beam  evaporation.  A  14-/xm- thick  layer  of  photoresist 
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FIG.  1.  (Color  online)  (a)  Cross-section  view  of  the  filter  board 
and  trap  chip  assembly,  (b)  Photograph  of  filter  board  and  mounted 
surface-electrode  trap,  (c)  Micrograph  of  the  central  region  of  the 
surface-electrode  trap.  The  electrodes  are  labeled.  The  direction  of 
the  external  quantization  field  |Bo|  ~  21.28  mT,  located  in  the  y-z 
plane,  is  shown.  The  ^^Mg+  ions  are  trapped  at  a  distance  J  ~  30  /xm 
above  the  surface,  at  the  position  labeled  trap  center. 

(Shipley  SPR220-7^)  is  spin  coated  onto  the  wafer,  exposed 
through  a  patterned  chromium  mask  by  contact  lithography, 
and  developed.  The  photoresist  pattern  defines  4.5-/xm-wide 
gaps  between  electrodes  in  the  central  region  of  the  trap;  the 
gap  width  is  increased  to  25  /xm  at  distances  :^2  mm  away 
from  the  trap  center.  Gold  is  then  electroplated,  in  a  140-mL 
sulfite  gold  solution  (Transene  TSG-250),  up  to  a  10.5-/xm 
thickness  onto  individual  approximately  22.5  x  22.5  mm^ 
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dies.  The  use  of  a  small  electroplating  volume  allows  us 
to  replace  the  solution  for  each  electroplating  process,  thus 
reducing  the  possibility  of  contamination.  After  removing  the 
photoresist  with  acetone,  the  seed  and  adhesion  layers  are  wet 
etched,  and  a  total  of  nine  (approximately  7.5  x  7.5  mm^)  trap 
chips  are  obtained  by  dicing.  To  provide  good  thermal  contact, 
the  chip  is  bonded  with  a  thin  layer  of  silver-filled  vacuum 
compatible  epoxy  (EPO-TEK  H21D)  to  a  solid  copper  support 
[Eig.  1(a)].  This  support  structure  also  holds  a  printed  circuit 
board  for  in- vacuum  filtering  of  the  control  potentials.  We 
use  low-pass  RC  filters  with  R  =  1  k^2  (ANAREN,  part  no.: 
H2B15081001E2LO)  and  C  =  820  pE  (NOVACAP,  part  no.: 
0504N821 JIOIP).  In  addition,  the  printed  circuit  board  is  used 
for  connecting  both  rf  electrodes  to  a  resonant  quarter-wave 
step-up  transformer  (gj-es  —  350,  when  loaded  with  the  trap) 
[46].  Connections  between  all  trap  chip  electrodes  and  the 
filter  board  are  made  with  gap-welded  gold  ribbons  with 
a  50  X  500  /xm^  cross  section  and  length  ^2  mm.  Each 
microwave  electrode  is  connected  to  a  microstrip  line  on  the 
filter  board,  which  is  soldered  to  a  SMA  jack  on  the  input  end. 
The  other  end  of  each  microwave  electrode  is  shorted  to  ground 
at  the  edge  of  the  chip  via  gap- welded  gold  ribbons.  In- vacuum 
coaxial  cables  connect  the  vacuum  feedthroughs  with  the 
SMA  connectors.  All  of  the  assembly  steps  are  performed 
in  a  cleanroom  environment  to  minimize  the  presence  of 
contamination,  e.g.,  dust  particles  on  the  surface  of  the  trap 
chip,  which  can  generate  unwanted  conduction  paths  and  stray 
electric  fields,  possibly  enhanced  when  exposed  to  UV  laser 
beams. 

An  rf  peak  voltage  Vrf  15  V  to  60  V  at  ^2rf  — 
In  X  71.6  MHz  is  applied  to  both  rf  electrodes,  providing  the 
x-z  radial  confinement  of  the  ions  at  a  distance  J  30  /xm 
above  the  surface.  The  rf  frequency  is  stabilized  to  the 
center  of  the  quarter-wave  resonance  of  the  resonator  with  a 
feedback  loop  to  the  frequency  modulation  input  port  of  the 
rf  generator.  The  six  control  electrodes  (Cl  to  C6)  are  biased 
between  —10  V  and  -h  10  V  to  provide  confinement  along  the 
trap  axis.  The  potentials  are  applied  by  six  digital-to-analog 
converters  (DACs;  for  details  see  Sec.  III).  Single-ion  motional 
frequencies  range  between  /axM  =  0.5  MHz  and  2.0  MHz 
in  the  axial  direction,  while  two  nondegenerate  frequencies 
between  /radial  =  3  MHz  and  12  MHz  aligned  along  adjustable 
directions  (Sec.  IV  C)  in  the  radial  x-z  plane  are  used. 

We  calculate  the  trapping  potential  /trap  in  the  radial  and 
axial  directions  by  means  of  the  gapless  plane  approximation 
[45];  /trap  is  generated  by  combining  the  rf  pseudopotential 
and  the  control  potential. 

With  Vrf  =  35  V  and  Vci-C6  =  {-0.801,0.641,-0.801, 
0.750,-0.384,0.750}  V,  Eig.  2(a)  shows  the  calculated 
/trap(-^.^)  in  the  radial  plane  around  the  trap  center  at  x  =  y  = 
^  =  0,  while  Eig.  2(b)  illustrates  /trap(T)  in  axial  direction. 
The  overall  trap  depth,  limited  by  a  saddle  point  in  the  x-z 
plane  [Eig.  2(a)],  is  :^13  meV  for  this  choice  of  potentials. 
We  denote  the  two  radial  center-of-mass  modes  by  their 
relative  frequencies,  low-frequency  (LE)  and  high-frequency 
(HE)  mode;  their  orientations  are  indicated  by  arrows  in 
Eig.  2(a).  The  orientation  of  these  modes  is  determined  by  the 
control  potentials  while  the  rf  pseudopotential  is  cylindrically 
symmetric  in  a  small  volume  around  the  trap  center.  Due  to 
the  tapered  geometry  of  the  rf  electrodes,  we  find  a  weak 
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EIG.  2.  (Color  online)  Calculated  trap  potential  /^p  (in  units 
of  meV)  as  a  function  of  distance  from  the  trap  center  (rf  null 
and  axial  minimum),  which  is  located  30  /xm  above  the  surface. 
Applied  potentials  are  given  in  the  text,  (a)  Potential  in  radial  plane 
at  y  =  0.  The  dot  marks  the  trap  center.  The  vectors  indicate  the 
calculated  motional  mode  orientation.  Here  q^lf  denotes  the  angle 
between  the  LE-radial  mode  vector  and  the  x  axis,  (b)  Potential  in 
the  axial  direction  at  x  =  z  =  0.  The  contribution  from  the  axial  rf 
pseudopotential  (see  text)  is  scaled  by  a  factor  of  10  for  visibility. 


axial  rf  pseudopotential,  with  a  curvature  corresponding  to 
a  motional  frequency  /axial  —  20  kHz.  The  applied  axial 
confinement,  however,  is  dominated  by  the  control  potentials 
and  its  frequency  is  simulated  to  be  /axM  —  0.51  MHz. 
Additionally,  we  notice  the  presence  of  significant  stray 
potentials,  which  must  be  taken  into  account  in  the  simulations 
to  reproduce  the  observed  motional-mode  directions  and 
frequencies  (Sec.  IV  C). 

We  employ  a  set  of  “shim”  fields,  designed  through  simula¬ 
tions  to  produce  orthogonal  displacements  of  the  ion(s)  along 
X,  y,  and  z  while  approximately  preserving  the  control-field 
curvatures,  for  micromotion  nulling  (Sec.  IV  B)  and  to  obtain 
oscillating  magnetic-field  maps  (Sec.  IV  A).  Due  to  the  limited 
DAC  resolution  (Sec.  Ill),  the  ion  position  can  be  controlled  to 
a  precision  of  {Ax,Ay,Az}  —  {16,37,5}  nm  for  typical  trap 
frequencies  /radial  —  7.0  MHz  and  /axial  —  1.5  MHz. 

A  static  magnetic  field  of  magnitude  |Bo|  ~  21.3  mT  is 
aligned  parallel  to  the  trap  surface  at  an  angle  of  —15°  with 
respect  to  the  z  axis  [Eig.  1(c)]  and  provides  the  ions’  internal- 
state  quantization  axis.  At  this  field  strength,  the  ground 
state  ^Si/2  \F  =  3,mF  =  1)  =  ||)  to  \F  =  =  1)  =  It) 

hyperfine  transition  frequency  /o  —  1.686  GHz  is  first-order 
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EIG.  3.  (Color  online)  Ground- state  hyperfine  structure  of 
(nuclear  spin  5/2)  at  |Bo|  —  21.3  mT.  The  microwave  driven 
hyperfine  transitions  used  to  prepare,  detect,  and  control  the  qubit 
states  are  depicted  as  arrows.  The  cooling  and  detection  laser  beam 
is  indicated  as  a  (blue)  wavy  arrow.  (Inset)  Illustration  of  the  qubit 
states  (bold  lines)  and  their  radial  motional  states  (gray).  The  “red” 
and  “blue”  sideband  transitions  as  well  as  the  qubit  transition  are 
indicated  with  arrows. 

insensitive  to  magnetic  field  changes,  while  the  quadratic 
frequency  deviation  is  :^233  kHz/(mT)^;  we  use  the  |t)  and 
II)  states  as  the  qubit.  Here,  F  is  the  total  angular  momentum 
and  m/7  is  the  projection  of  the  angular  momentum  along 
the  magnetic  field  axis.  Figure  3  illustrates  the  ground-state 
hyperfine  structure,  highlighting  the  qubit  states.  We  track 
I  Bo  I -field  drifts  every  20  to  30  min  with  a  Ramsey  experiment 
on  the  1 3 , 3 )  -to- 1 2, 2)  transition  (excited  with  microwave  fields), 
which  has  a  magnetic-field  sensitivity  of  :^19.7  MHz/mT 
and  adjust  a  computer-controlled  correction  coil  current 
accordingly.  We  observe  day-to-day  |Bo| -field  variations  of 
less  than  0.02  mT. 

Two  overlapping  a + -polarized  laser  beams  for  Doppler 
cooling  on  the  ^  5’i/2  1 3 , 3 )  ->  /2 1 4, 4)  cycling  transition  prop¬ 
agate  parallel  to  the  magnetic  field.  A  first  beam  (BDD),  ap¬ 
proximately  —240  MHz  detuned  from  this  transition,  performs 
initial  Doppler  cooling  and  optical  pumping.  A  second  beam 
(BD)  is  red  detuned  from  this  transition  by  approximately  half 
the  natural  line  width  (r/2  ~  27r  x  20  MHz)  and  is  used  for 
final  Doppler  cooling  and  state  preparation  into  the  1 3, 3)  level 
of  the  ^5’i/2  ground  state. ^  State  preparation  and  cooling  takes 
1  ms;  during  the  final  50  /xs,  only  the  BD  beam  is  applied. 
For  state  detection,  the  BD  beam  is  tuned  to  resonance  for 
approximately  200  /xs,  discriminating  the  |3,3)  state  from  the 
other  hyperfine  ground  states  through  resonance  fluorescence. 
The  fluorescence  photons  are  detected  by  a  photon-multiplier 
tube  (PMT)  detector.  In  most  experiments,  the  1 3, 3)  population 
is  transferred  via  two  microwave  preparation  pulses  to  ||) 


^Note  that  here  we  choose  the  |3,3)  state  for  preparation  instead  of 
the  |3,— 3)  state  used  in  Ref.  [32]. 


to  initialize  the  qubit.  To  detect  the  ||)  state,  two  pulses 
are  used  to  transfer  the  ||)  population  to  the  |3,3)  (bright) 
state,  while  the  |t)  population  is  “shelved”  into  the  |2,  — 1) 
(dark)  state  (Fig.  3).  Shelving  suppresses  repumping  of  the  1 1 ) 
state  into  the  bright  state  during  detection  and  enhances  the 
contrast  between  |t)  and  ||)  by  approximately  10%.  During 
the  detection  period,  the  PMT  registers  counts  for  the 
bright  state  and  ~0.2  counts  for  the  dark  state  on  average  per 
ion.  Stray  light  accounts  for  ~0. 1  counts  on  average.  The  BDD 
and  BD  beams  are  both  focused  into  the  vacuum  chamber  with 
a  10-cm  focal  length  lens,  and  their  powers  are  adjusted  to  be 
1  and  ~30  /xW,  respectively;  both  have  waists^  of  15  /xm  at 
the  position  of  the  ions.  To  create  ions,  a  third  continuous- wave 
(CW)  laser  beam  (denoted  PI:  A.  285  nm,  power  :^3  mW, 
and  beam  waist  ~  15  /xm)  is  directed  antiparallel  to  the  cooling 
beams  and  is  used  to  photoionize  neutral  Mg  atoms  inside 
the  trapping  region.  An  oven  tube,  holding  enriched  ^^Mg,  is 
mounted  approximately  2  cm  away  from  the  trap  center.  After 
allowing  the  oven  to  heat  up  for  30  s,  the  PI  is  applied  for  a  few 
seconds  until  ion(s)  are  loaded.  The  oven  current  is  adjusted  to 
a  loading  rate  of  less  than  one  ion  per  second,  so  the  PI  beam 
can  be  turned  off  manually  after  the  desired  number  of  ions 
are  detected  on  a  CCD  camera.  At  a  background  pressure  of 
2  X  10“^  Pa,  the  storage  period  of  a  single  ion  is  observed  to  be 
up  to  a  few  hours,  with  Doppler  cooling  applied  periodically 
every  1  or  2  ms.  Typical  trapping  lifetimes  for  two  ions  are  on 
the  order  of  30  min. 

III.  MICROWAVE  SYSTEM 

To  implement  microwave  control  of  the  ion(s)  spin  and 
motional  states,  we  built  a  system  assembled  mostly  from  com¬ 
mercially  available  discrete  microwave  components  (Fig.  4). 
For  hyperfine  state  manipulation  and  single-qubit  operations, 
magnetic  fields  from  microwave  currents  in  electrode  MW2 
drive  transitions  at  frequencies  between  1.2  to  2.3  GHz, 
spanning  all  hyperfine  transitions  in  the  ^5*1/2  manifold.  As 
noted  above,  we  use  a  series  of  two  transfer  pulses  (carrier 
TT  pulses)  to  transfer  population  from  |3,3)  to  |3,1)  via  |2,2) 
to  initialize  the  qubit  state.  We  have  demonstrated  tt  times  for 
qubit  spin  flips  as  short  as  :^20  ns  [32] ;  however,  here  we  adjust 
the  power  of  the  microwave  signals  to  levels  corresponding  to 
TT  times  on  the  order  of  a  few  microseconds,  limited  by  the 
timing  resolution  of  our  data  acquisition  system  (see  below). 

To  drive  “blue”  (BSB)  and  “red”  (RSB)  sidebands  of  the 
qubit  transition  (Fig.  3),  we  apply  currents  simultaneously 
to  MWl,  MW2,  and  MW3  at  /s  =  /o  +  /radial  or  /s  =  /o  - 
/radial-  The  amplitudes  and  phases  of  the  three  signals  are 
adjusted  to  minimize  the  oscillating  magnetic  field  at  the 
position  of  the  ion(s),  while  maximizing  the  field  gradient 
(Sec.  IV  A).  This  control  is  enabled  by  the  six  phase  shifters 
and  six  attenuators. 

The  microwave  setup  is  integrated  with  a  data  acquisition 
system,  controlled  by  a  field-programmable  gate  array  (FPGA, 
Xilinx  Virtex  IV)  with  a  timing  resolution  of  16  ns  [48].  The 
entire  system  comprises  16  digital  (TTL)  output  channels, 
24  DAC  channels,  eight  direct-digital  synthesizer  (DDS) 


^Denotes  the  half- width. 
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EIG.  4.  Schematic  diagram  of  the  microwave  system  used  to 
control  the  microwave  fields  applied  to  the  ion(s).  It  enables  individual 
control  of  phases  and  amplitudes  of  the  currents  in  all  three  microwave 
electrodes  for  two  frequencies  =  /q  =b  /radial  (BSB  and  RSB) 
around  1.69  GHz,  within  a  range  of  about  200  MHz.  It  also  enables 
a  signal  between  1.2  and  2.3  GHz  applied  to  MW2  for  hyperfine 
state  control  and  single-qubit  operations.  The  gray  components  are 
computer  controlled. 

modules,  and  two  digital  (TTL)  inputs  for  registering  PMT 
counts.  Each  DDS  (Analog  Devices  AD9858)  can  generate 
signals  at  frequencies  to  400  MHz  with  a  frequency  resolution 
of  0.233  Hz  and  phase  resolution  of  0.022°.  The  DACs 
(National  Instruments  NI  PXI-6733)  have  an  update  rate  of 
500  kHz  and  can  be  set  between  —10  V  and  +10  V  with 
a  resolution  of  0.305  mV.  We  use  six  of  these  channels  to 
apply  potentials  to  the  trap  control  electrodes  and  14  channels 
to  control  microwave  components.  In  addition,  a  fast,  two- 
channel  DAC  is  used  to  generate  arbitrary  waveforms  with  a 
50-MHz  update  rate,  a  voltage  range  from  — 10  V  to  10  V,  and 
a  resolution  of  0.305  mV.  This  DAC  is  programed  via  USB 
and  triggered  by  the  data  acquisition  FPGA. 

We  use  three  DDS  modules  as  sources  for  three  frequency 
octupling  modules,  each  consisting  of  three  doubling  stages 
(Mini-Circuits,  two  FD-2+  and  one  FK-3000+)  in  series  with 
commercial  amplifiers  and  filters.  Two  of  these  generate  the 
RSB  and  BSB  signals,  which  pass  through  individual  pulse 
shaping  and  control  stages  (phase  and  amplitude  control). 
The  pulse  envelopes  for  the  BSB  and  RSB  signals  are 
produced  with  a  multiplier  (Analog  Devices  ADL5391)  used 
to  shape  pulse  amplitudes  and  controlled  by  one  channel  of 
the  fast  DAC.  In  addition,  a  voltage-controlled  attenuator 
(Mini-Circuits  ZX73-2500)  adjusts  the  overall  power  of  a 
pulse.  After  this  pulse-shaping  stage,  the  signals  are  split  and 
directed  to  three  individual  voltage-controlled  phase  shifters 
(:^12°/V,  custom  made  by  Pulsar  Microwave  Corporation) 
and  attenuators  (:^1.6  dB/V,  Mini-Circuits  ZX73-2500)  to 
control  the  power  and  phase  of  the  signals  delivered  to  the 
individual  microwave  electrodes.  Before  the  RSB  and  BSB 
signals  are  combined  onto  the  inputs  of  each  microwave 
electrode,  they  are  amplified  by  high-power  amplifiers  (Mini- 
Circuits  ZHF-30W-252-S+).  Following  three  high-power 
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combiners  (Werlatone  D7630),  three  coaxial  cables  connect 
the  microwave  setup  with  the  vacuum  feedthroughs.  At  the 
inputs  to  the  feedthroughs,  we  measured  the  frequency  noise 
of  the  signals  at  1.6865  GHz  and  estimated  a  negligible  effect 
on  the  observed  fidelity.  The  third  octupling  module  generates 
frequencies  to  drive  the  hyperfine  carrier  transitions,  and  its 
signal  is  combined  with  one  of  the  sideband  lines  attached  to 
electrode  MW2.  Its  power  is  adjusted  with  fixed  attenuators  to 
yield  carrier  tt -pulse  durations  of  a  few  microseconds. 

On  a  daily  basis,  we  calibrate  the  frequencies  and  tt -pulse 
durations  for  the  transfer  and  single-qubit  pulses.  From  day 
to  day,  we  observe  fractional  variations  of  ^10“^  for  the 
pulse  durations.  Due  to  our  quantization  field  calibrations 
and  adjustments,  the  hyperfine-transition  frequency  values 
stay  constant  within  a  statistical  uncertainty  of  approximately 
0.3  kHz. 

IV.  EXPERIMENTAL  PROCEDURES  AND 
CALIBRATION  SEQUENCES 

We  use  a  single  ion  to  calibrate  experiments  and  set 
parameters  for  the  two-qubit  gate.  After  an  ion  is  loaded  into 
the  trap,  micromotion  along  the  projection  of  the  detection 
laser  beam  in  the  radial  plane  (z  direction)  is  minimized  by 
measuring  the  variation  in  laser  fluorescence  from  the  BD 
beam  [47]  as  a  function  of  the  “horizontal”  -shim  field  (Fig.  5). 
Since  the  trap  rf  quadrupole  is  tilted  in  the  radial  plane  by 
aRF  —  —20°  from  the  x  axis  (Fig.  5),  the  horizontal-shim 
field  is  chosen  to  push  along  a  vector  rotated  by  2aRF 
from  the  z  axis.  A  push  along  this  axis  will  result  in  a 
maximal  variation  of  the  horizontal  micromotion  amplitude. 
We  verify  this  compensation  field  with  a  nulling  scheme  using 
microwave  magnetic-field  gradients  (described  in  Sec.  IV  B). 
This  technique  can  also  be  used  to  minimize  the  micromotion 
in  the  vertical  (y)  direction.  We  find  that  the  trapping  fields 


rf  quadrupole  electric  field 


z(pm) 

FIG.  5.  (Color  online)  Calculated  trap  rf  quadrupole  electric  field 
in  the  radial  plane.  The  quadrupole  axes  are  rotated  by  q^rf  ~  —20°. 
The  arrows  indicate  the  horizontal  (red)  and  the  vertical  (gray) 
micromotion  component.  Only  the  horizontal  component  can  be 
detected  by  the  detection  laser  beam.  The  red  and  gray  lines  indicate 
the  two  shim  directions  used  for  the  compensation  of  horizontal  and 
vertical  micromotion. 
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are  perturbed  by  stray  fields,  with  corresponding  potential 
0stray  The  exact  origin  of  these  fields  is  unknown,  but  they 
are  consistent  with  a  contact  potential  patch  and  can  be 
compensated  for  as  described  in  Appendix  A. 

Following  micromotion  nulling,  we  calibrate  the  mi¬ 
crowave  magnetic  fields  used  for  driving  the  radial  sidebands 
of  the  qubit  transition  (Sec.  IV  A).  From  these  experiments, 
we  extract  the  orientation  and  the  strength  of  the  microwave 
magnetic-field  gradient.  As  a  last  calibration  step  to  set  up  the 
two-qubit  gate  experiment,  we  align  the  radial-mode  vectors 
relative  to  the  microwave  magnetic  field  by  rotating  the  modes 
with  an  additional  control  field  (Sec.  IV  C). 


A.  Controlling  and  mapping  the  microwave  near  field 

To  implement  motional  sideband  transitions  at  frequencies 
/s,  a  gradient  of  the  microwave  magnetic  field  is  necessary. 
However,  any  oscillating  field  at  the  ion(s)  position  will  cause 
off-resonant  carrier  transitions  and  ac  Zeeman  shifts  [21], 
both  of  which  will  inhibit  precise  control.  Thus,  the  magnetic 
field  amplitude  should  be  suppressed  as  much  as  possible  at 
the  position  of  the  ion.  For  this  purpose,  there  is  a  single 
combination  of  relative  current  amplitudes  and  phases  in 
electrodes  MW  1 ,  MW2,  and  MW3  that  provides  an  oscillating 
magnetic  field  gradient  at  the  ions  without  an  oscillating 
magnetic  field.  In  principle,  we  can  determine  this  combination 
through  simulation,  but  we  realize  this  experimentally  by 
adjusting  the  phases  and  amplitudes  of  the  three  microwave 
currents  to  minimize  ac  Zeeman  shifts  imposed  on  the  ion.  The 
experimental  sequence  is  shown  in  Fig.  6(a). 

The  nulling  procedure  relies  on  the  measurement  of 
ac  Zeeman  shifts  on  superpositions  of  selected  hyperfine 
states  \a)  and  \b).  In  a  coordinate  system  with  unit  vectors 

=  Bo/ 1  Bo  I,  and  055=02x0^,  the  magnetic-field 

components  oscillating  at  /§  can  be  written  as 
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where  B\\  and  are  the  strengths  of  the  oscillating  magnetic- 
field  components  parallel  and  perpendicular  to  the  quantization 
field  Bo.  Here  ^  is  the  rotation  angle  of  the  perpendicular 
magnetic-field  component  polarization  ellipse  with  the  x  axis 
and  y  is  a  phase.  For  6  =  0,  the  oscillating  magnetic  field  in 
the  plane  perpendicular  to  Bo  is  linearly  polarized.  For  6  = 
ib7r/4,  the  perpendicular  field  component  has  pure  circular 
polarization.  The  B\\  field  couples  levels  with  Am/?  =  0, 
whereas  the  B^  field  drives  transitions  with  Am/?  =  zb  1 .  These 
fields  induce  an  ac  Zeeman  shift  5/acz  =  c^B\-\-  c\\  on  the 
\a)  ^  \b)  transition,  where  the  coefficients  may  depend  on  fs 
and  6.  Note  that  we  implicitly  include  the  Bloch-Siegert  shifts 
(energy  shift  resulting  from  the  counter-rotating  part  of  the 
field)  in  the  calculation  of  the  coefficients  in  5/acz- 

First,  consider  \a)  =  |2,0)  and  \b)  =  |3,0).  A  B±^ 
component  applied  near  /o  is  red  detuned  from  the 
|2,0)^{|3,-1),|3,1)}  and  |3,0)  ^  {|2,-1),  |2,1)}  tran¬ 
sitions,  thereby  increasing  the  energy  of  |a)  =  |2,0) 
and  decreasing  the  energy  of  jZ?)  =  |3,0).  From  the 
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FIG.  6.  (Color  online)  Realization  of  optimal  microwave-field 
configurations  for  driving  motional  sideband  transitions,  (a)  Pulse 
sequence  to  detect  the  presence  of  residual  oscillating  magnetic 
fields  at  the  ion(s)  position  through  phase  accumulation  from  the 
corresponding  ac  Zeeman  shift,  while  the  ion  is  in  a  superposition 
state  of  |a)  and  \b).  The  spin  echo  n  pulse  makes  the  result  insensitive 
to  magnetic -field  fluctuations  slower  than  the  time  scale  of  an 
individual  experiment  (typically  ~500  /xs).  We  observe  minimal 
phase  pickup  when  the  oscillating  magnetic  held  is  minimized  at 
the  position  of  the  ion.  (b)  Extension  to  measure  the  ac  Zeeman 
shift  as  a  function  of  displacement  Sr  (see  text),  (c)  Experimental 
data  mapping  out  the  ac  Zeeman  shift  in  the  radial  x-z  plane  by  use 
of  the  superposition  state  of  \a)  =  |2,0)  and  \b)  =  |3,1).  The  fltted 
isolines  are  obtained  from  the  quadrupole  oscillating  magnetic  held 
component  B\\  along  Bq  through  its  resulting  ac  Zeeman  shift  Vacz, 
which  is  proportional  to  5^.  The  angle  q^mw  is  the  rotation  angle  of 
the  microwave  quadrupole  held,  but  the  directions  of  minimum  and 
maximum  values  of  |55||  /Sr  \  are  along  angles  2q'mw  and  2q'mw  +  90° 
[Fig.  8(a)]. 


Clebsch-Gordan  coefficients  for  these  transitions,  we  de¬ 
termine  c±_  (0.357  —  0.092  sin  26)  Hz/(/xT)^  at  /§  =  /o.  A 

B\\  component  near  /o  is  red  detuned  from  the  |3,0)  ^ 

1 2,0)  transition  and  also  increases  the  \a)  ^  \b)  transition 
frequency:  cy  0.468  Hz/(/xT)^  at  /§  =  /q.  Both  and  cy 
are  positive,  of  similar  magnitude  for  all  6  and  depend  only 
weakly  on  /§  around  /o  within  the  useful  range  of  motional 
sideband  frequencies,  so  this  choice  of  \a)  and  \b)  can  be  used 
to  minimize  both  B±  and  B\\ . 

Another  useful  choice  is  \a)  =  |2,0)  and  \b)  =  |3,1). 
Here,  the  dominant  ac  Zeeman  shift  results  from  B\\ 
coupling  to  the  |3,1)  ^  |2,1)  transition:  cy  ^  —49.08 
/[(  fs  -  /o)/MHz]  Hz/(AtT)2,  and  cx  (0. 199  + 
0.220  sin  2e)Hz/(/iT)^  at  = /q.  This  transition  can 
be  used  to  minimize  B\\  even  further  due  to  its  large  cy  value 
by  use  of  currents  in  electrodes  MW2  and  MW3  (those 
electrodes  contribute  dominantly  to  B\\  due  to  their  geometry). 
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We  produce  the  data  of  Fig.  6(c)  using  this  transition  and  an  ex¬ 
perimental  sequence  presented  in  Fig.  6(b).  With  Bmw  =  0  at 
position  8r  =  =  {0,0,0}  in  the  coordinate  system 

of  Fig.  1(c)  [not  the  frame  of  Eq.  (1)]  we  expect  to  first  order 

Bmw  =  B'  cos(27tfy) 


^  cos(2q:mw) 

0 

sin(2o'Mw)  ^ 

/8x\ 

0 

0 

0 

Sy 

y^sin(2Q'Mw) 

0 

—  cos(2q'mw)  j 

\Sz) 

Here  B'  and  q^mw  characterize  the  strength  and  orientation 
of  the  quadrupole  for  small  displacements  8r  <  1.5  /xm  from 
the  field  null  point.  By  fitting  the  ac  Zeeman  shift  resulting 
from  the  projection  B\\  of  Bmw  onto  Bq  to  the  experimental 
data  of  Fig.  6(c),  we  extract  B'  =  35.3(4)  T/m  and  q'mw  = 
—26.5(7)°.  From  this  fit,  we  determine  a  small  shift  of  the 
center  of  the  magnetic  quadrupole  with  respect  to  the  trap 
center  by  |5x,5z}  {—0.11,0.23}  /xm,  limited  by  the  level  of 

control  over  the  relative  phases  and  amplitudes  in  MW  1 ,  MW2, 
and  MW3.  We  detect  phase  accumulations  over  a  time  = 
250  /xs;  we  cannot  measure  phase  shifts  smaller  than  tt/IO 
during  that  time.  This  limits  the  ac  Zeeman  shift  sensitivity  to 
~0.2  kHz  and  leads  to  finite  fields  at  the  trap  center.  This  is  also 
the  level  of  control  of  the  experimental  day-to-day  variations 
for  the  nulled  configuration  are  ^0.016  dB  for  the  relative 
powers  of  MW2  and  MW3  with  respect  to  MWl  and  ^0.12° 
for  the  relative  phases  of  MW  1  and  MW 3  with  respect  to  MW2. 

Numerical  simulations  of  the  microwave  currents  in  the 
three  electrodes  by  use  of  ANSYS  HFSS  software  [46]  provide 
a  value  of  ^simiA  — ^4  T/m  for  a  current  amplitude  of 
1.0  A  in  conductor  MWl,  while  the  amplitudes  and  phases 
of  currents  in  MW2  and  MW3  are  adjusted  to  give  a  nulled 
field  configuration  at  the  position  of  the  trap  rf  null.  We  find 
the  following  current  amplitude  ratios  R  and  phase  shifts 
AO  in  MW2  and  MW3  relative  to  the  current  in  MWl  for 
the  nulled  configuration:  {/^mw2,^mw3}  —  {0.82,0.96}  and 
{AOmw2, AOmws}  —  {161.6°,  — 14.6°}.  A  comparison  with 
the  experimental  current  ratios  is  difficult  because  of  slightly 
varying  attenuations  in  the  signal  paths  leading  to  the  three 
microwave  electrodes.  However,  from  simulations,  we  extract 
an  angle  —25°,  close  to  the  value  fitted  to  the  exper¬ 

imental  data  (also  Sec.  IV  B).  A  similar  level  of  agreement 
between  simulation  and  experiment  is  found  in  Ref.  [40]. 
Furthermore,  Ref.  [40]  discusses  an  optimized  trap  chip  design 
based  on  numerical  simulations.  In  our  experiments,  current 
densities  are  simulated  to  be  as  high  as  :^2  x  10^^  A/no?  near 
the  edges  of  the  conductor  MW  1 ,  with  a  total  current  of  :^0.8  A 
in  MWl,  for  the  above  stated  experimental  settings,  producing 
a  gradient  of  35  T/m.  We  estimate  the  corresponding  power 
dissipation  in  all  electrodes  due  to  all  currents  from  the  bulk 
resistance  of  gold  (—2.2  x  10“^ ^2  m)  in  the  central  region  of 
the  chip  (±200  /xm  along  y)  to  be  ~80  mW.  In  Fig.  7  the  current 
density  distribution  in  electrode  MW  1  and  the  induced-current 
densities  in  the  neighboring  electrodes  RFl  and  C2  are  shown. 
The  simulated  magnetic  and  electric  microwave  fields  for  the 
nulled  configuration  are  plotted  in  Fig.  8.  The  influence  of  the 
microwave  electric  fields  on  the  trapping  potential  is  discussed 
in  Appendix  B. 
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EIG.  7.  (Color  online)  Cross-section  view  of  the  simulated 
microwave  current  density  distribution  in  the  trap  electrodes  when  the 
magnetic  field  is  nulled  at  the  ion  position  (only  showing  electrodes 
REl,  MWl,  and  C2).  The  densities  shown  here  correspond  to  a 
gradient  of  ~  35  T/m.  Currents  in  microwave  electrodes  induce 
currents  in  neighboring  electrodes,  here  illustrated  for  REl  and  C2. 

B.  Microwave  approach  for  rf  micromotion  nulling 

In  the  presence  of  excess  rf  micromotion,  the  microwave 
magnetic  near-field  gradient  can  also  drive  micromotion- 
induced  sidebands  on  carrier  transitions  [47] .  This  effect  can  be 
used  to  null  micromotion  along  any  direction  in  the  radial  x-z 
plane  [refer  to  the  frame  of  Fig.  1(b)].  Therefore,  it  may  have 
an  advantage  over  certain  optical  schemes  that  are  sensitive  to 
micromotion  only  along  the  direction  of  a  laser  beam  (or  along 
the  direction  of  the  effective  k  vector  for  stimulated-Raman 
transitions)  [47].  For  example,  in  a  surface-electrode  trap,  the 
direction  perpendicular  to  the  surface  can  be  difficult  to  probe 
with  laser  beams,  because  scattered  light  from  the  trap  surface 
may  interfere  with  detection  and/or  cause  electrode  charging. 

The  micromotion  can  be  written  as  (ymm^x  ±  yMM^^  ± 
ZuM^z)  cos(27r/RF0.  In  the  presence  of  a  microwave  magnetic 
field  gradient  of  frequency  /§  and  in  the  coordinate  frame 
of  Fig.  1(c),  the  ion  in  its  rest  frame  experiences  fields  with 
frequencies  /§  ±  /rf  and  amplitude  Bion,  cf.  Eq.  (2): 


^  cos(2q:mw) 

0 

sin(2Q:Mw)  \ 

^xmm  ^ 

0 

0 

0 

0 

sin(2Q'Mw) 

0 

—  cos(2q'mw)  / 

\zmm  / 

(3) 

For  /s  =  /hfs  3=  /rf,  where  /hfs  is  the  frequency  of  a 
transition  between  states  in  the  hyperfine  manifold,  the 
ion  experiences  an  oscillating  field  at  frequency  /hfs-  The 
resulting  Rabi  rate  is  proportional  to  |  Bion  I  and,  therefore, 
to  tmm,  the  micromotion  amplitude.  To  minimize  the  excess 
micromotion,  we  move  the  ion  in  the  radial  plane  (along 
the  shim  directions  illustrated  in  Fig.  5)  to  a  position  that 
minimizes  the  Rabi  rate  of  the  micromotion  sidebands.  The 
oscillating  field  Bion  has  projections  B\\  and  B±  on  Bq. 
Therefore,  we  can  sense  micromotion  along  any  direction  in 
the  x-z  plane,  since  it  is  always  detectable  by  either  sl  Am f  =0 
or  Amp  =  ±1  transition  (or  both).  In  the  experiment,  we 
use  the  |2,0)  ^  |3,1)  and  |3,1)  ^  |2,1)  transitions,  with  a 
magnetic-dipole  matrix  element  :^0.414  /xr  and  1.001  /xr, 
respectively,  where  /xr  denotes  the  Bohr  magneton. 

To  estimate  the  sensitivity  of  this  nulling  method,  we 
assume  B'  =  35  T/m  and  a  transition  matrix  element  for 
/hfs  of  /^B-  To  further  simplify  the  estimate,  the  rotation 
angles  of  the  microwave-magnetic  quadrupole  [auw  —  —30°, 
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EIG.  8.  (Color  online)  Illustration  of  the  simulated  microwave 
fields  in  the  radial  plane,  (a)  Magnetic  quadrupole  field  as  a  function 
of  displacement  from  the  trap  center,  (b)  Electric  field.  Due  to 
the  pseudopotential  from  the  electric  field  the  ion(s)  are  displaced 
from  the  trap  rf  null  as  indicated  in  gray.  The  pseudopotential  from 
the  electric  field  displaces  the  ion(s)  from  the  trap  rf  null  (see 
Appendix  B).  The  estimated  ~0.4  nm  displacement  for  typical  trap 
parameters  is  shown  in  gray,  scaled  by  a  factor  10^  for  visibility. 


Fig.  8(a)]  and  rf-electric  quadrupole  («rf  —  —20°,  Fig.  5) 
are  assumed  to  be  oriented  such  that  in  the  rest  frame  of 
the  ion,  we  achieve  a  maximum  oscillating  field  component 
that  drives  the  transition  (q^mw  =  Ck^RF/2).  For  the  given 
misalignment,  this  approximation  is  good  to  within  15%  for 
the  estimated  sensitivity.  Note  that  in  the  experiments  both 
components  of  the  microwave-field  gradient  can  be  used,  and 
the  misalignment  does  not  limit  the  sensitivity  of  this  method. 
The  Rabi  rate  for  the  micromotion  sideband-induced  carrier 
transition  is  then  given  by  [21] 

1  ^  /Xg 

^MM  =7^^  ^MM  •  (4) 

2  2n 

For  our  parameters,  and  assuming  a  tt/IO  rotation  detection 
sensitivity  on  the  micromotion  sideband  during  a  drive  time 
Tp  =  2  ms,  a  micromotion  amplitude  of  tmm  =  0.1  nm  can  be 
detected. 

In  addition  to  detecting  micromotion,  we  use  this  method 
to  probe  characteristics  of  the  microwave-magnetic  and 
rf-electric  field  by  measuring  the  micromotion  sideband  Rabi 
rate  as  a  function  of  position.  Figure  9  shows  the  Rabi 
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FIG.  9.  Rabi  rate  on  the  red  micromotion  sideband  of  the  1 2,0)  ^ 
1 3, 1)  transition  as  a  function  of  ion  position.  The  isolines  result  from 
a  model  fit  to  the  data,  from  which  5'  =  35.1(4)T/m  and  q^mw  = 
—31.1(2)°  are  extracted  from  Eq.  (3). 

rate  on  the  red  micromotion  sideband  of  the  |2,0)  ^  |3,1) 
transition  versus  ion  position.  Using  these  data  together  with 
the  calculated  shape  of  the  rf  quadrupole,  we  extract  the 
parameters  B'  and  q'mw  of  the  microwave  quadrupole  field, 
finding  B'  =  35.1(4)T/m  and  Ck^MW  =  —31.1(2)°.  This  map 
was  taken  with  a  different  realization  of  the  trap  described 
in  Sec.  IV  A  but  of  the  same  design,  so  imperfections  in 
construction  and  small  variations  in  the  electrode  impedance 
may  have  affected  the  measurement.  Differences  in  impedance 
affect  the  strength  of  the  gradient,  so  the  agreement  in  B'  values 
is  probably  fortuitous.  However,  we  find  again  a  value  for  q'mw 
that  is  close  to  the  simulated  value  ~  —25°  (Sec.  IV  A). 

C.  Control  of  the  orientation  of  the  radial  motional  modes 

After  determining  the  characteristics  of  the  microwave 
quadrupole,  we  align  the  orientation  of  the  radial  motional 
modes  to  the  microwave  field  gradient  to  optimize  sideband 
Rabi  rates.  Since  the  qubit  is  controlled  with  a  Amp  =  0  tran¬ 
sition,  its  motional  sidebands  are  driven  by  8B\\/8r  [Fig.  8(a)]. 
To  overlap  the  HF-mode  vector  with  55||/5r  Uax,  we  use 
an  additional  control  potential  0rot  with  potentials  V(5f.c6  = 
{-1.694,2.298,-1.694,0.311,0.501,0.311}  V  applied  to  the 
control  electrodes,  which  leads  to  a  mode  orientation  q^lf  —  5° 
in  the  absence  of  any  additional  control  potential.  In  the 
experiments,  the  total  potential  0totai  =  0trap  +  0stray  +  ^rot0rot 
breaks  the  symmetry  of  the  rf  pseudopotential  and  sets  the 
orientation  of  the  radial  modes.  The  value  of  the  scaling 
coefficient  Vj-ot  is  used  to  adjust  the  orientation  of  the  radial 
modes.  For  the  given  total  potential  0totai  including  the  stray 
potential  (Appendix  A),  negative  Vj-ot  rotates  the  LF-mode 
vector  toward  q^lf  ^  0°  and  for  positive  values  it  rotates 
toward  q^lf  ^  —90°.  Furthermore,  0i-ot  is  designed  such  that 
the  ion  position  remains  unchanged  while  all  trap  frequencies 
change  as  a  function  of  the  scaling.  A  comparison  of  the 
simulated  frequency  changes  with  the  actual  experimental  data 
can  be  used  to  determine  and  adjust  the  orientation  of  the  radial 
modes.  In  Fig.  10(a)  such  a  comparison  between  simulated 
and  measured  mode  frequencies  shows  good  agreement. 
Figure  10(b)  illustrates  the  corresponding  mode-angle  values 
extracted  from  the  simulation.  A  similar  comparison  between 
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EIG.  10.  (Color  online)  Control  over  the  radial  mode  orientation, 
(a)  Comparison  of  simulation  (lines)  with  experimental  data  of 
radial-mode  frequencies  (error  bars  are  too  small  to  be  visible)  as 
a  function  of  the  scaling  of  the  control  field,  (b)  Mode  angle  q^lf 
extracted  from  the  simulation.  At  an  angle  q^lf  =  ^q^mw  —  —53° 
the  high-frequency  mode  is  aligned  along  and  couples 

maximally  to  the  microwave-field  gradient. 

simulation  and  experimental  data  is  performed  for  different 
trap  geometries  in  Ref.  [49] .  Note  that  for  the  simulation  we 
had  to  include  0stray  to  find  agreement  with  the  experimental 
data,  details  are  given  in  Appendix  A.  When  choosing 
Ck^LF  =  2aMw  —  —53°,  we  observe  the  maximal  coupling  of 
the  microwave-field  gradient  to  the  HE  mode,  while  the  LF 
mode  couples  very  weakly  to  the  microwave  field.  We  find  a 
corresponding  Rabi  rate  for  driving  the  HF  sideband  transition 
of  the  motional  ground  state  of  —In  x  2.0  kHz,  with  the 
microwave-field  gradient  adjusted  to  35  T/m. 

V.  PRESENT  LIMITATIONS  EOR  TWO-QUBIT  GATES 

Internal-to-motional  state  coupling  as  well  as  an  entangling 
two-qubit  gate  based  on  the  simultaneous  application  of 
detuned  red  and  blue  sideband  drives  [42]  was  presented 
in  [32].  Using  this  gate,  entangled  states  with  a  fidelity  of 
approximately  0.76  were  produced.  In  these  experiments  we 
used  an  out-of-phase  radial  mode  of  two  ions,  because  it  is 
less  sensitive  to  electric-field  noise  than  the  corresponding 
center-of-mass  mode.  The  out-of-phase  motional  heating  was 
measured  to  be  0. 2-0.5  quanta/ms,  while  the  center-of-mass 
motional  heating  rate  was  approximately  5  quanta/ms  (for  ra¬ 
dial  mode  frequencies  around  7  MHz).  Heating  corresponding 
to  an  increase  of  a  single  quantum  of  motion  on  the  mode 
used  for  an  entangling-gate  operation  effectively  decoheres 
the  gate  operation  [50].  Since  in  [32]  the  entangling  gate 
durations  were  approximately  250  /xs,  even  the  out-of-phase 
mode  heating  can  cause  substantial  gate  errors.  To  investigate 
this  effect  for  our  experimental  conditions,  we  numerically 
solve  the  master  equation  given  in  Ref.  [50].  We  determine 
that  this  effect  alone  limits  the  fidelity  to  approximately  0.8. 
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Furthermore,  we  observed  drifts  of  the  motional  sideband 
frequencies  in  the  experiments,  which  may  be  caused  by  mode 
frequency  drifts  and/or  varying  ac  Zeeman  shifts  from  drifts 
of  the  residual  magnetic  and  electric  microwave  fields.  We 
estimate  the  frequency  variation  during  the  gate  pulse  to  be 
less  than  0.5  kHz  from  these  effects,  while  a  typical  detuning 
for  the  gate  is  around  4.5  kHz.  This  systematic  effect  can 
contribute  an  additional  5%  loss  of  fidelity. 

VI.  CONCLUSION 

In  summary,  we  have  given  details  of  the  surface-electrode 
trap  and  microwave  electronics  used  for  the  microwave  near¬ 
field  quantum  control  experiments  described  in  Ref.  [32] .  We 
have  described  the  experimental  procedures  used  to  control, 
shape,  and  calibrate  the  microwave  fields  that  implement 
motional- sideband  transitions  using  the  near-field  gradients. 
We  have  introduced  a  novel  technique  for  rf  micromotion 
detection.  Furthermore,  we  presented  and  experimentally 
verified  a  method  to  rotate  the  radial  mode  axes.  Finally, 
we  investigated  limitations  in  fidelity  of  the  two-qubit  gates 
in  Ref.  [32]  through  simulations  and  examined  the  impact 
of  leading  error  sources.  We  determined  that  the  motional 
heating  during  the  gate  pulse  is  the  dominant  source  of 
infidelity.  A  significant  reduction  (by  a  factor  of  100  or  more) 
of  the  heating  rate  should  lead  to  entangled  state  fidelities 
higher  than  0.9.  A  reduction  of  ambient  heating  may  be 
achieved  by  surface  treatments  of  the  electrode  structure  as 
described  in  Refs.  [38,39]  and/or  cryogenic  cooling  of  the 
electrodes  [35-37].  In  addition,  higher  microwave  currents 
can  increase  the  gate  speed  and  thereby  can  reduce  the  effect 
of  motional  heating.  The  corresponding  increased  heat  load 
may  be  mitigated  by  reducing  electrode  resistance  through 
low-temperature  operation. 
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APPENDIX  A:  COMPARISON  OE  TRAP  SIMULATION 
WITH  MEASUREMENTS 

To  compare  the  trapping  potential  calculations  from  Sec.  II 
to  the  experiment,  we  first  determine  the  electric  field  strengths 
at  the  position  of  the  ion  to  compensate  the  rf  micromotion 
(Sec.  IV).  We  then  measure  the  mode  frequencies  and 
find  the  orientation  of  the  radial  mode  vectors  (Sec.  IV  C), 
characterized  by  a\jp  (Fig.  2).  We  find  a  stray  field  Estray  — 
720  X  {-0.54,0.10,-0.84}  V/m  and  a  stray  field  curvature 


013437-9 


U.  WARRING  et  al 


PHYSICAL  REVIEW  A  87,  013437  (2013) 


tensor,  given  by  the  Hessian  matrix  of  0stray  at  the  center  of  the 
trap, 
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After  including  0stray  into  the  potential  simulations,  we  find 
reasonably  good  agreement  (motional  frequencies  and  mode 
angles  agree  to  within  10%  or  better)  between  the  simulations 
and  the  measurements  (Sec.  IV  C)  for  several  different  settings 
of  rf  voltage  and  control  potentials.  The  stray  potential  might 
be  caused  by  work  function  variations  on  the  electrode 
surface  [5 1-54]  (e.g.,  from  varying  Mg  contamination)  or  stray 
charges  localized  in  the  central  region  of  the  trap  chip.  One  way 
of  modeling  the  stray  effects  is  a  rectangular  patch  with  size 
Ay  =  35  /xm  and  Az  =  5  /xm  located  on  the  chip  surface, 
whose  center  is  displaced  by  z  =  —25  /xm,  y  =  +3  /xm 
relative  to  the  trap  center  and  that  is  biased  to  —  1 . 15  V  (typical 
work  function  differences  can  exceed  one  volt).  We  observe 
temporal  variations  in  the  stray  potential  effects,  perhaps  due  to 
uncontrolled  charging  of  the  electrodes  which  may  be  caused 
by  laser  light  impinging  on  the  electrodes  and/or  changing 
Mg  coverage  of  the  chip  surface  during  loading  of  the  trap. 
In  particular,  the  applied  compensation  fields  vary  within  a 
few  tens  of  volts  per  meter  on  a  day-to-day  basis,  and  the 
axial-mode  frequency  can  vary  more  than  200  kHz  over  the 
course  of  a  few  weeks  with  many  loading  cycles. 

We  estimate  the  deviation  of  ion  displacements  from 
simulated  values  due  to  stray-electric  fields  to  be  less  than 
a  few  percentages  for  typical  displacements  from  the  rf  null 
of  zbl  /xm  at  radial-mode  frequencies  of  :^7.0  MHz.  We  do 
not  include  this  effect  in  our  analysis  of  the  experiments 
described  in  Sec.  IV  A  and  Sec.  IV  B.  However,  for  the 
calibrations  described  in  Sec.  IV  C,  we  take  the  stray-electric 
field  into  account  because  of  its  significant  impact  on  the  radial 
motional-mode  vector  directions. 


APPENDIX  B:  INFLUENCE  OF  THE  MICROWAVE 
ELECTRIC  FIELDS  ON  THE  TRAPPING  POTENTIAL 

In  addition  to  the  oscillating  magnetic  field  Bmw,  an 
electric  field  Emw  (at  the  microwave  frequency)  is  also 
present  [Fig.  8(b)].  If  the  null  points  for  Emw  and  Bmw  fields 
differ,  the  pseudopotential  produced  by  Emw  can  push  the 


ions  away  from  the  main  trap  rf  pseudopotential  null.  From 
numerical  simulations,  we  determine  that  the  ion  position  in 
a  7.0  MHz  trap  is  shifted  by  {5v,5z}  —  {0.3,— 0.2}  nm  when 
B'  ~  35  T/m  [Fig.  8(b)].  This  displacement  is  much  smaller 
than  the  precision  for  positioning  the  ion  with  the  shim  fields 
(Sec.  II)  and,  therefore,  this  shift  cannot  be  resolved. 

Another  systematic  effect  is  caused  by  the  curvature  of  the 
microwave  pseudopotential.  It  modifies  the  mode  frequency 
(based  on  simulations)  by  about  0.5  kHz,  depending  on  the 
mode  direction.  This  can  be  a  significant  shift  for  the  two-qubit 
gate.  However,  when  calibrating  the  sideband  frequencies  this 
effect  is  included,  since  we  probe  the  sidebands  with  the 
microwave  drive,  which  inherently  includes  the  microwave 
pseudopotential  produced  by  the  gate  pulse.  To  experimentally 
check  for  this  shift,  we  measure  the  mode  frequencies  by 
applying  a  oscillating  potential  (coherent  excitation)  to  one 
of  the  control  electrodes,  while  the  ion  is  in  the  |3,3)  state. 
We  then  detect  the  ion  fluorescence  as  a  function  of  excitation 
frequency:  a  decrease  in  the  resulting  fluorescence  indicates 
that  a  mode  of  the  ion’s  motion  has  been  resonantly  excited 
[55].  We  compare  the  measured  mode  frequencies  with  and 
without  the  microwave  fields  while  applying  the  excitation 
pulse.  With  the  microwaves  applied,  we  observe  shifts  of  the 
radial-mode  frequencies  of  :^3  kHz.  The  observed  shift  is 
about  7  times  higher  than  what  we  expect  from  our  estimate; 
this  discrepancy  is  not  understood.  The  mode  frequency  shift 
(for  a  7.0-MHz  trap)  for  the  ion  in  the  |3,3)  state  due  to  the 
spatially  varying  ac  Zeeman  shift  is  estimated  to  be  less  than 
:^10  mHz  and  can  be  neglected  here. 

We  checked  for  motional  excitations  when  simultaneously 
applying  two  tones  near  /o  zb  /radiah  as  required  for  the  two- 
qubit  gate  (Sec.  V).  For  pulse  lengths  up  to  1  ms  (4  times  longer 
than  typical  gate  times)  we  observe  no  change  in  the  motional 
state  within  an  uncertainty  of  approximately  zbO.5  quanta. 
However,  when  we  apply  the  two  tones  at  /o  zb  /radial /2,  we 
record  coherent  motional  excitations  corresponding  to  a  few 
tens  of  quanta  after  50  /xs.  The  origin  of  the  observed  excitation 
may  be  the  gradients  of  the  two  microwave  pseudopotentials 
(with  difference  frequency  /radial)-  We  estimate  this  effect  to 
be  negligible  for  our  experimental  conditions.  However,  we 
cannot  rule  out  that  the  excitation  is  caused  directly  by  an 
electric  field  at  frequency  /radial-  This  field  may  originate  from 
mixing  of  the  two  microwave  tones  due  to  nonlinear  behavior 
of  one  or  more  of  the  microwave  elements. 
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